High-order synchronization, transitions, and competition among 
Arnold tongues in a rotator under harmonic forcing 
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We consider a rotator whose equation of motion for the angle 6 consists of the zeroth and first 
Fourier modes. Numerical analysis based on the trailing of saddle-node bifurcations is used to 
locate the n:l Arnold tongues where synchronization occurs. Several of them are wide enough for 
high-order synchronization to be seen in passive observations. By sweeping the system parameters 
within a certain range, we find that the stronger the dependence of 6 on 8, the wider the regions 
of synchronization. Use of a synchronization index reveals a vast number of very narrow n:m 
Arnold tongues. A competition phenomenon among the tongues is observed, in that they "push" 
and "squeeze" one another: as some tongues widen, others narrow. Two mechanisms for transitions 
between different n:m synchronization states are considered: slow variation of the driving frequency, 
and the influence of low-frequency noise on the rotator. 

PACS numbers: 05.45.Xt, 45.20.dc, 05.40.Ca, 87.19.Hh 
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I. INTRODUCTION 

When two or more oscillatory processes are coupled, 
there exists the possibility of their becoming synchro- 
nized. Where their autonomous frequencies are different 
but close, synchronization is understood as the adjust- 
ment of those frequencies as a result of coupling. Even 
when such systems operate on different timescales, syn- 
chronization may still appear as an adjustment of their 
frequencies to an integer ratio, an effect known as high- 
order synchronization or synchronization of order n:m. 

Synchronization of order n:m has been extensively 
studied, both experimentally and theoretically (cf. 1| 
for a review). For instance, Simonet and co-workers 2\ 
investigated a ruby nuclear magnetic resonance laser with 
delayed feedback. The undriven laser exhibited periodic 
oscillations of light intensity at a frequency vq ~ 40 Hz. 
An external periodic voltage, either sinusoidal or square- 
wave, was then added in the feedback loop. In both cases, 
synchronizations of different order n:m were observed. 
Another example is the electrical rotator consisting of a 
Josephson junction, shunted by a capacitor, and fed with 
a constant external current. Experiments show that this 
system can be synchronized when a periodic driving cur- 
rent is applied, or where two junctions are coupled [3|. 
There are many additional examples. 

Synchronization has also been observed extensively in 
biology. One example is the cardiorespiratory system, 
considered in the pioneering works by Kenner et al, Hildc- 
brandt and Raschke Schafer and co-workers Q pro- 
posed the synchrogram as a tool to visualize cardiorespi- 
ratory synchronization. When plotting the instantaneous 
respiratory phase at the occurrence of a heartbeat ver- 
sus time, they found horizontally striped plots for some 
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subjects, thereby revealing n:l synchronization between 
heart and respiration. Toledo et al [6| showed that the 
probability of such synchronization happening by chance 
was extremely small. In measurements on anaesthetised 
rats, Stefanovska et al Q observed lengthy synchroniza- 
tion epochs, and transitions from one ratio to another. 
They suggested that such transitions might be useful in 
monitoring depth of anaesthesia. 

There are several nonlinear models yielding mm syn- 
chronization. Arnold proposed Q a map of the cir- 
cle into itself, obtaining synchronization "tongues" and 
calculating their widths in the approximation of small 
coupling. In fact, any orientation-preserving homcomor- 
phism h-.S 1 — > S 1 of the circle into itself presents such 
regions of n:m locking [nj. Integrate-and-fire models 
provide several examples of Arnold tongues [ll|. High- 
order synchronization regions were also obtained by Glass 
and Sun [12| for an impulse-driven Poincare oscillator. 
Schilder and Peckham [13| treated Arnold tongues nu- 
merically, and they obtained tongues for the system of 
two coupled Van der Pol oscillators; here the tongues are 
quite narrow, so that the probability of locking in a real, 
noisy system, would therefore be very small. Simonet 
and co-workers proposed a model that reproduced the 
synchronizations observed in the laser. For cardiores- 
piratory synchronization, Kotani et al 1 141 developed a 
model, based on those of DcBoer et al [15[ and Seidel 
and Herzel [l(J, which is supported on both physiologi- 
cal and mathematical principles. The model involves a 
somewhat complicated system consisting of several oscil- 
lators and interactions, and includes some non-analytic 
parts (integrate-and-fire). Because of technical difficul- 
ties encountered when tackling nonlinear models, the 
deep mechanisms through which synchronization takes 
place are not yet understood in general, so that there 
is still no way of predicting which equations, and which 
values of their parameters, will or will not yield synchro- 
nization. 
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The present paper has two main purposes. First we 
present a systematic study of high-order synchronization 
in a particularly simple system, a rotator under harmonic 
forcing, for which we can establish the roles played by 
each of its parameters in synchronization. Although we 
study a specific system, and although the deep mecha- 
nisms responsible for synchronization are not unveiled in 
this paper, we report below two results that we believe 
will be useful in the quest for those deep mechanisms: 
we show that the main Arnold tongues are wider when 9 
depends more strongly on 0; and that competition occurs 
between synchronization regions. Our second purpose is 
to discuss and explore two possible mechanisms giving 
rise to an extensively observed phenomenon: transitions 
between the different n:m synchronization ratios. The 
mechanisms considered here arise from time variability. 
We will also discuss briefly why time variability hinders 
the analysis of synchronization in experimental data. 

The paper is organized as follows: In Sec. |TT] we in- 
troduce the simple rotator whose synchronization prop- 
erties are to be considered. Sec. [HI] discusses how its re- 
gions of n:l and n:m synchronization are identified and 
reports the main results obtained under stationary condi- 
tions, including the observation of competition between 
the tongues. Time variability and its effect on transitions 
between different synchronization states is discussed in 
Sec. IIVI Sec. IVl summarizes the main conclusions. 



II. THE SYSTEM 

The generic equation for the angle 9 of a rotator with- 
out external interaction is 

9 = f(0), 

where / is a 27r-periodic function [l7|. Therefore, such 
a system can be studied systematically by considering 
functions / up to a certain number k of harmonics, and 
allowing a bigger number k + 1 of harmonics at the next 
stage of the study. In this paper, we start the study for a 
function consisting of the zero harmonic (the "constant 
force" , thanks to which a rotator has the features of a self- 
sustained oscillator) and the first harmonic. By means of 
a translation in the value of 9, it can always be written 
as 



= ag + a± cos 9, 



(1) 



with ao and a\ constants, where (JXJ) is an Adler-type 
equation [l9l |. We assume ao > |ai| so that, in the ab- 
sence of interaction (nothing is added to IJTJ), the angle 
9 continuously increases. For oo much bigger than oi, 
9 increases at an almost constant rate. For a\ close to 
ao, however, 9 varies strongly with 9. We show that (fT]) 
can synchronize to an external forcing, exhibiting a wide 
variety of Arnold tongues, and we discuss the processes 
that may be responsible for transitions between different 
n:m synchronization ratios. 



The equations for the overdamped pendulum, and the 
overdamped Josephson junction, are of just this type. 
The equation of motion of a pendulum driven by a con- 
stant torque K is described by 



, K 
* + 7* + n sin* = — , 



(2) 



where k is the frequency of small oscillations, 7 > is 
the damping constant, and / is the moment of inertia. 

In a resistively shunted Josephson junction, the current 
is a sum of three contributions: a superconducting cur- 
rent I c sin\I/, a current V/R through the resistance, and 
a capacitance current V C . The parameter I c is called 
the critical supercurrcnt of the junction. The relation- 
ship between the potential V and the "angle" ^ is given 
by the Josephson formula 



(3) 



where e is the electronic charge. Summation of the three 
terms yields the equation for a junction fed with an ex- 
ternal current I 



I — I c sin * 



h d* Chd 2 ^ 



2eR dt 2e dt 2 



(4) 



Equations @ and @ coincide. In the overdamped limit, 
when the term with the second derivative can be ne- 
glected - in the case of the Josephson junction, this 
means that there is no capacitor in the circuit - these 
equations reduce to the one that we study in this paper 
((T|) . Thus the results that we obtain below will be appli- 
cable to the overdamped pendulum and the overdamped 
Josephson junction. 

As Q is analytically integrable, we find that the fre- 
quency for the non-interacting rotator is 



(5) 



We now consider the effect of an external harmonic 
force on the rotator 



= ao + ci\ cos 9 + B sin(w t), 



(6) 



with B > 0. This equation applies to a number of situa- 
tions in nature. In the case of a Josephson junction, Eq. 
([6]) describes when the system is fed with a continuous 
intensity plus a harmonic one. In relation to cardiores- 
piratory synchronization, |6]) can be regarded as a very 
simple model in which the rotator (TTJ) models the heart, 
and the addition of the harmonic component models its 
interaction with respiration [2(| • Eq. © has been studied 
[l| for large amplitudes B of the harmonic component. 
However, we are interested here in the Arnold tongues 
down to very small driving amplitude, in order that we 
can also apply the results to the weakly coupled rotator. 

Eq. ^ defines a circle map with the following prescrip- 
tion: let us call U the time at which the external driving 
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is at its i th maximum. We define 0, = 0(ij). As ([6]) is 
a first-order differential equation and we have an initial 
condition 0(t,), we could integrate to obtain 0(f). Let 
ti + \ be the time at which the external driving is at its 
i + 1 maximum, and let 0j + i = 0(ij + i). That is how we 
have the map h-.S 1 — > 5 1 defined as 0;+i = h(6i). We 
would therefore expect the existence of Arnold tongues 

Bill- 



ill. SYNCHRONIZATION AND ARNOLD 
TONGUES 

A. Synchronization of a rotator driven by an 
external periodic force 

Now we introduce the concept of n:m synchronization 
that will be used in this work. Suppose we have a rotator 
driven by a Tg-periodic external action 

9 = f{9) + g(t), with 

/(0 + 2tt) = /(0) V0, g(t + T E )=g(t)Vt. 

Then, any T-periodic motion of 0, 6(t+T) = 9(t) mod 2tt, 
must have a period that is a multiple of the driving pe- 
riod. Let to G N be such that T = mT E . Let n G NU{0} 
be the number of times that the angle crosses = mod 
2tt with > in one of its periods T = mTg. We then 
say that the rotator is n:m synchronized to the external 
driving g. The synchrogram consists then of horizonal 
lines, albeit not necessarily equally spaced. Note that 
any periodic motion of the forced rotator automatically 
implies n:m locking to the external action for some n and 

TO. 



B. Saddle- node bifurcation 

Putting the above definition of mm synchronization 
into mathematical terms, we say that the rotator gets 
synchronised n:m to the external force if there exists a 
stable root (zero) for the function 

h nm (9)=h m (9)-2%n-6, (7) 

where h m stands for the return map h composed with 
itself m times. In general, if h nm has two or more roots, 
there is at least a stable fixed point. Regardless of the ini- 
tial condition, the trajectory is attracted towards a stable 
fixed point. The function h nm depends, of course, on the 
parameters of ©, so it changes as we vary the external 
driving frequency u>. At the moment of transition from 
two zeros to no zero, there is only one root, on which 
h nm is tangent to the horizontal axis, this single zero is 
a half-stable fixed point. At this moment, a saddle-node 
bifurcation takes place. 

Obtaining the borders of the Arnold tongues therefore 
involves retrieving the two driving frequencies at which 



the saddle-node bifurcation takes place for different val- 
ues of the driving amplitude B. For this purpose a con- 
tinuation software was written in C. The first step was to 
obtain the Poincare return map h, so the interval [0, 2tt] 
was divided into many points. At each of these, the func- 
tion h m was obtained by integration of the differential 
equation ^ with the fourth-order Runge-Kutta method, 
thereby obtaining h nm (9) = h m (9) — 2irn — 9. For given 
values of B and u, we know that we are inside the tongue 
if h nm has two or more roots, and outside if h nm has no 
zeros. Hence, for a given value of B, we can trace the 
left and right values of u> at which the saddle-node bi- 
furcation takes place, up to the desired precision. We 
start with zero driving amplitude, B = 0; for this value 
the tongue consists of only one point u = to uq/u - this 
will not be forced when we study the bifurcations in the 
"flexible" way for very narrow tongues, see Section [ill Fl 
For the next value of the coupling (positive but close to 
0), we start from u> — m luq/h and look for the two bi- 
furcation points. Then, for sequentially increasing values 
of B, we first guess approximately the bifurcation points 
on the left and right boundaries by linear extrapolation 
from the two former bifurcation frequencies in each case. 
Starting from this guessed value, the programme looks 
for the correct bifurcation point. In some tongues and 
for small values of B, the synchronization region is so 
narrow that the programme cannot retrieve the bifur- 
cation points. In such cases, the programme skips this 
value of B, goes to the next B, and sets the starting point 
(the "guessed value" of the bifurcation frequency) based 
on the former bifurcation points that were successfully 
retrieved. 



C. The n:l synchronization regions 



We obtained the first n:l tongues for three combina- 
tions of parameters (ao,ai) that, according to (|5j). yield 
the same autonomous frequency for the rotator i/q ~ 
1.0811: (7.08116, 2), (9.85714, 7.14286) and (16.46637, 
15). Note that dimensionless units are used for all vari- 
ables throughout this paper. The results are plotted in 
Fig. [TJ The tongues with to > 1 are very narrow and the 
programme fails to retrieve them by the method of seek- 
ing the saddle-node bifurcations. Later we will use other 
less accurate methods for the high-order tongues. As of- 
ten reported in the literature Q, the widest tongues are 
those for m = 1. For our system ©, as n increases, the 
widths of the to = 1 tongues decrease and the tongues 
become closer to one another. A point of interest here 
is that there are several Arnold tongues that are wide 
enough to be seen, not only by seeking them in a tuned 
experiment, but also through passive observations in the 
presence of noise, e.g. the cardiovascular system [j| 0,113] . 

In our second parameter set, a\ is bigger than in the 
first set, and longer again in the third set compared to 
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FIG. 1: (colour online) Arnold tongues for the externally- 
driven rotator (JBJ), with parameters (ao,ai) equal to: 
(7.08116, 2) in (a); (9.85714, 7.14286) in (b); and (16.46637, 
15) in (c). The frequency v and amplitude B of the external 
driving are plotted on the abscissa and ordinate axes respec- 
tively. The tongues are ml, with n increasing from 1 to 9 as 
we move from right to left. 



though the frequency v$ remains the same. The effect on 
the Arnold tongues is immediately evident on comparing 
the three plots in Fig. [TJ the tongues become markedly 
wider, thereby favouring synchronization. We may there- 
fore conjecture that a strong dependence of the instan- 
taneous frequency on angle (for the uncoupled rotator) 
may help high order synchronization. 

In order to check this conjecture more thoroughly, we 
fix the value of the driving amplitude B to 1.5, and we 
compute the widths of the tongues at B = 1.5 for dif- 
ferent values of the parameters ao and a\ . We swept ao 
from 2 to 30, and a\ from 2 to (almost) ao. In order to 
compute the width of the tongues, the code was modified, 
so that the two driving frequencies at which the saddle- 
node bifurcations take place were computed only for this 
value of B. Starting from u> = muio/n, we trailed the 
bifurcation points first with small precision, then with a 
greater one, recursively up to the desired precision. We 
plot in Fig. [5] the widths of the tongues at B = 1.5 versus 
the autonomous frequency vq and the parameter a\ (to 
which ao is related through In the figure we see 

that, for the same autonomous frequency the widths 
of the tongues dramatically increase with oi, confirming 
the conjecture that a strong dependence of 8 on 8 (for 
the uncoupled rotator) favours synchronization. Regard- 
ing this conjecture, it would be very interesting to add 
higher harmonics to (fTJ) and to study how the widths of 
the tongues for the driven rotator change. Nevertheless, 
in this case, the frequency vq of the isolated rotator is 
not an elementary function, so it would perhaps be more 
convenient to calculate this frequency numerically. 



D. Use of the synchronization index 

As mentioned above, the tongues with m > 1 are very 
narrow and the programme fails to retrieve them by the 
method of seeking the saddle-node bifurcations. In order 
to obtain the synchronization regions in such cases, we 
have been obliged to use other (albeit less accurate) ap- 
proaches. The first one is via a synchronization index: we 
generated data by numerical integration of Eq. ^ with 
the fourth-order Runge-Kutta method; and we analysed 
them as though they were experimental data, by com- 
puting their synchronization index. With the definition 
of synchronization we mentioned at the beginning, it is 
natural to quantify it by how close the following index 
11411 is to zero 



Hm.n. — 



1 



2tt m (L — n) 



E 

i— n-fl 



27rm) 



(8) 



the second one. Thus, the value of 8 depends much more 
strongly on 8 as we move from (a) to (c) in Fig. [TJ al- 



where L is the total number of crossings of 8 = mod 
2-7T with 8 > 0, and <& is the phase of the external action 
when the i th crossing occurs. 
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FIG. 2: (color online) Widths of the Arnold tongues as functions of parameters a± and uq: (a) 1:1; (b) 1:2; (c) 1:3; and (d) 1:4. 
They are expressed in units of the driving frequency u, at driving amplitude B — 1.5. Note the differences in ordinate scale. 



E. The n:m synchronization regions 

To illustrate the n:m synchronization regions, we again 
plot the driving frequency v on the horizontal axis, and 
the coupling or driving amplitude B on the vertical axis. 
The horizontal axis was divided into 5,000 points, and 
the vertical one into 30 points. For each of the 150,000 
corresponding pairs of points, the Eq. © was integrated. 
The first 60 zero-positive crossings were discarded to re- 
move transients. The maximum accepted value of the 
index H mn ([5]) for a point to be considered as belonging 
to the n:m Arnold tongue was 0.0005. 

Fig.[3]plots the results obtained for two of the parame- 
ter sets. We only plotted those tongues whose horizontal 
widths were larger than 0.00026 for at least one value 
of the coupling B (an arbitrary cutoff chosen by trial 
and error, in order to plot neither too many nor too few 
tongues) [23]. We can see in that figure the vast variety 
of synchronization regions that the system yields. There 
are even some on this Fig. with m = 10, the highest value 



of m that we trailed. As we reduce the lower cutoff of the 
width criterion, many new tongues arise in the plot. The 
effect of synchronization with the external action can be 
inferred from the fact that these tongues are not vertical 
lines, but curved. 

Furthermore, comparison of Figs. [22a) and(3Jb) shows 
that the tongues become closer to vertical straight lines 
for smaller ax, implying that synchronization is more and 
more a matter of fine tuning of the external driving fre- 
quency rather than of the interaction of the rotator with 
the external action. This phenomenon provides further 
confirmation of our hypothesis that a strong dependence 
of instantaneous frequency on angle (for the uncoupled 
rotator) helps high-order synchronization. 

F. Competition among the tongues 

The curving of the high-order tongues in Fig. [3]Ja) can 
be attributed to the occurrence of competition among 
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FIG. 3: (colour online) Detail of the synchronization regions 
inside a small interval. The parameters (ao,ai) are: (a) 
(9.85714, 7.14286); (b) (7.08116, 2). The tongues are, from 
left to right: (a): 5:1, 33:8, 41:10, 4:1, 39:10, 35:9, 31:8, 27:7, 
23:6, 19:5, and 34:9; (b): 5:1, 29:7, 33:8, 37:9, 41:10, 4:1, 
39:10, 35:9, 31:8, and 27:7 



tongues: in this case the 4:1 tongue widens as the cou- 
pling increases and "pushes" the high-order tongues. Our 
circle map h is always invertible, because equation © 
can be integrated with reversed time. Therefore, Lemma 
1.5 in fl(| says that the Arnold tongues cannot overlap. 
That is the basis on which we argue that there must be 
competition among the tongues. 



We note here several signs of such competition. First, 
the 37:9 tongue appears in Fig. [3^b) but not in [3ja), 
because it is narrower there. We can interpret this effect 
as 33:8 and 41:10 squeezing it for the set of parameters 
used in[3Ja). 

In the same way, the widths for 33:8 and 41:10 obtained 
via the synchronization index in Fig. EJa) are of the order 
of 10~ 4 at B -> 0, and of the order of 10~ 5 at B = 3. 
Nevertheless, the width of the synchronization region as 
B — > must tend exactly to zero. As indicated above, 
the computation of Arnold tongues via a synchronization 
index (using (JSj) in our case) is less accurate. Because of 
this evident contradiction between the obtained width 
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FIG. 4: (colour online) (a) Tongues 4:1 (right) and 33:8 
(second from the right) for parameters (ao, ai)=(9. 85714, 
7.14286). For clarity, the width of the 33:8 tongue in (a) 
is shown as larger than it really is: this plot is just intended 
to show the location, (b) Width of the 33:8 tongue, plotted 
as a function of the driving amplitude B. The location of 
the 33:8 tongue in (a), and its width in (b), were obtained by 
trailing the bifurcations in their "flexible" meaning, using a 
threshold of 0.001. 



of the tongue as B — > (width bigger than zero) and 
the actual value of the width (strictly zero), we need a 
more defined approach. The first try was to obtain the 
bifurcation points with the first programme but it was 
impossible: although we set a high precission, we failed 
to find any zeros for the function h nm . 

That is why we are introducing the third method of 
obtaining approximately Arnold tongues for tricky cases 
like this one: trailing the bifurcations in their "flexible" 
meaning. Instead of considering that we are inside the 
tongue when h nm has a stable zero, we relax the condition 
and just ask h nm to have an "almost (half stable) zero" , 
defined as a point where h nm takes an absolute value 
smaller than a threhold (set here to be 0.001). It means 
that, if we measure the angle of the rotator at the times 
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when the external action has a maximum, this angle will 
stay for a long time around this minimum (or rather, min- 
imum in absolute value). What happens afterwards de- 
pends on the behaviour of h nm . If h nm increases quickly 
(in absolute value) when we separate from the minimum, 
the angle will slip and come back promptly to the min- 
imum (modulus 2tt). If h nm is still very small when we 
separate from the minimum, the angle will have a slow 
drift when iterating h nm so, in practice, the system will 
be regarded as synchronized, unless the measurements 
are very precise. 

The tongue 33:8 obtained this way is plotted in Fig. 
HJa), together with the 4:1 tongue obtained before. In 
Fig. EKb) we plotted the width of the 33:8 ton gue as 
a function of the driving amplitude B. Because of the 
"flexible" way of obtaining the tongue, its width tends 
to a value bigger than zero as B — > 0. The explanation 
is simple: without any driving (B = 0) there is no in- 
teraction, so the only possibility of having an apparent 
synchronization - not true synchronization, as there is 
no interaction - is to tune the external frequency to the 
value mvo/n. This frequency is the only point of the 
Arnold tongue in the limit B — > 0. The function h nm 
is there constantly equal to zero. If we make the exter- 
nal frequency slightly different from muo/n, with B = 0, 
the values returned by the function h nm will still be very 
small, so many of them will be bellow the threshold; that 
is the reason why the obtained width of the tongue in the 
limit B — ► is bigger than zero. Of course, this "flexi- 
ble" method is not completely accurate, and it gives the 
apparent contradiction above. Nevertheless, we are us- 
ing this method only to obtain qualitative conclusions, 
by comparing results all obtained with the same method. 
So we believe that the conclusions remain valid. 

As we can see from Fig. Hfb), the width of the 33:8 
tongue is more or less constant for B near zero, and then 
it starts dropping at around B — 1.5. Looking at Fig. 
IDJa), we may interpret B = 1.5 as the point where the 4:1 
tongue has approached 33:8 enough to start squeezing it. 
The width of 33:8 reaches a minimum at B = 2.5 and 
then it increases. In Fig. 0Ja), B = 2.5 is more or less 
the point where 4:1 switches from being concave (concave 
from outside) to convex, so it is consistent with the fact 
that the width of 33:8 starts increasing at this point: let 
us say that, before this point, 4:1 was "invading" 33:8; 
from this point, 33:8 "starts recovering from 4:l's inva- 
sion", so 33:8 gains space (width) and 4:1 switches from 
being concave to be convex (of course, there are, math- 
ematically speaking, an infinite number of tongues be- 
tween 33:8 and 4:1, but 4:1 is the most dominant in the 
area so it is the main influence). 

The cases discussed above reveal clear signs of the ef- 
fect of competition among tongues; the existence of such 
competition was predicted at the outset on strictly the- 
oretical grounds. 



IV. TIME VARIABILITY AND TRANSITIONS 
BETWEEN DIFFERENT SYNCHRONIZATION 
EPOCHS 

Transitions in time between states of different n:m syn- 
chronization ratios have often been observed experimen- 
tally MM- 

We now consider in turn two mecha- 
nisms that may give rise to such transitions in our system 
([5]): (i) variability of the driving frequency, and (ii) low- 
frequency noise. The latter is mathematically equivalent 
to slow variations of the autonomous frequency of the ro- 
tator. Both mechanisms correspond to time variability, 
and they can coexist. 

A. Variability of the driving frequency 

We again simulated the system ([5]), but modified such 
that we replaced Bsm(ujt) with B sin(<i>(£)), and varied 
the instantaneous frequency <fr(t) linearly in time 

$>(t) = at + b; $(f) = iai 2 + bt. (9) 

The choice of a linear variation was arbitrary, in the in- 
terests of clarity. With the parameters of ((6]) chosen to 
be (<z ,ai)=(9.85714, 7.14286), 5=4, we slowly swept 
$(£)/2tt from 0.1 to 1.3 in a total time of 72,000. Fig. 
[5]shows n:l synchronization epochs corresponding to the 
times when the external instantaneous frequency is in- 
side an Arnold tongue. As n decreases, the synchroniza- 
tion epochs last longer, and the transition regions widen. 
That is because, as mentioned above, the Arnold tongues 
are then wider and further separated. The n:2 epochs are 
also observed from the synchrograms: Fig. [6] shows 5:2 
and 3:2 synchronization epochs. 
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10" 3 Time 

FIG. 5: (colour online). An n:l synchrogram for the forced 
rotator while varying the driving frequency. Several n:l syn- 
chronization epochs can be observed. The straight line plots 
the instantaneous frequency of the external driving (right- 
hand ordinate axis). 
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FIG. 6: (colour online). Sections of an n:2 synchrogram for 
the forced rotator while varying the driving frequency, illus- 
trating a 5:2 synchronization epoch in (a) and a 3:2 epoch in 
(b). The straight line plots the instantaneous frequency of 
the external driving (right-hand ordinate axis). 

For a Josephson junction, transitions of this kind could 
be induced by variation of the frequency of the harmonic 
component of the current driving the system. Where we 
regard ^ as a simple model for the heart driven by the 
respiration, the transition mechanism would correspond 
to variation of the respiration rate. Although this rate 
can be consciously controlled [25| . people do not do so 
most of the time, and so the small natural variations in 
respiration rate can produce transitions in the sychro- 
nization ratio between it and the heart. 



B. Low- frequency noise 

The second transition mechanism considered here 
arises from the effect of the low-frequency noise that is 
often present in physical, biological and economic sys- 
tems. For instance, flicker noise [2(| occurs in almost 
all electronic devices. Another example is the cardiovas- 
cular system, where there are known to be metabolic, 



FIG. 7: (a) Transitions between 5:1 and 4:1 synchronization 
states, caused by the added low frequency component, (b) A 
9:2 epoch for time ~600 in the transition between 5:1 and 4:1. 

neurogenic and myogenic [13] low-frequency oscillatory 
processes. 

We therefore added to our system an extra low- 
frequency harmonic component, simulating Eq. ([6]) with 
(o ,ai)=(9.85714, 7.14286), B=4, w/2tt=0.213, with 
a term Lisin(wit) added to the RHS, Li=0.1 and 
cji/27t=0.001. This value of u>i may appear low compared 
to some natural systems, but the results we will obtain 
remain valid for larger w\ provided it is small compared 
with the frequencies of the other processes. Without 
adding the low frequency component, we would be inside 
the 5:1 tongue and the synchrogram would consist of 5 
perfect horizontal lines. As shown in Fig. \T[ a) , the corre- 
sponding simulation exhibits transitions between 5:1 and 
4:1. Note that we have added a harmonic low- frequency 
component; with low-frequency noise, the sequence of 
transitions will not be predictable in the manner seen 
here. In Fig. \T[b) we see a 9:2 epoch inside the transi- 
tion region between the two longest epochs. 

This second mechanism is mathematically equivalent 
to a slow variation of the intrinsic frequency of the rota- 
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tor. This can be seen by considering the low-frequency 
component to be absorbed inside a Q , yielding an ao, and 
therefore a vq, that slowly vary with time. In nature, 
however, the picture is completely different, as the low- 
frequency noise comes from outside the rotator. 

C. Time variability and synchronization 

Although we can still recognise synchronization epochs 
from the synchrograms by sight, time variability - the 
origin of the two transition mechanisms discussed above 
- means that synchrograms no longer consist of perfect 
horizontal lines. The test via a synchronization index, of- 
ten used when tackling real experimental data, must fail, 
even inside a small time window, unless we take a less 
strict threshold for the index. Our result is in agreement 
with that obtained in [22| by analysis of experimental 
data. In short, time variability interferes with synchro- 
nization; at least for our system and with the index ([5]) 
as a quantification of synchronization. A more detailed 
and quantitative study of the influence of time variability 
is in progress and will be reported elsewhere. 

V. CONCLUSIONS 

We have shown that the driven rotator Eq. ^ yields a 
vast variety of Arnold tongues. Our results are applica- 
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